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(N ! ABSTRACT 



We study a minimal one- loop radiative mechanism for generating small Majorana neutrino masses 
in inverse seesaw extensions of the Standard Model with two singlet fermions per family. The new 
feature of this radiative mechanism is that the one-loop induced left-handed neutrino mass matrix 
is directly proportional to the Majorana mass matrix of the right-handed neutrinos. This is a very 
economical scenario without necessitating the existence of non-standard scalar or gauge fields. 
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I. INTRODUCTION 



43. 

The understanding of the extra-ordinary smallness of neutrino masses, along with their observed large mixings [l[ 
O-l remains one of the major puzzles in the Standard Model (SM), and so a potential portal to New Physics. The simplest 
theoretical scenario which may explain this neutrino puzzle would be to introduce singlet Majorana masses that break 
the global (B — L)-symmetry of the SM. Within the SM, this can be parameterized through the non-renormalizable 
dimension-5 operator due to Weinberg [2j]: A.y (Lj&) (LJ<&) /A, where Li = (u i7 li)~[ (with i = e, /i,r) is the SU(2)l 
lepton doublet, $ = (0 + ,0°) T is the SM Higgs doublet, and A is an effective mass scale of the New Physics. After 
I ■ electroweak symmetry breaking via the vacuum expectation value (VEV) of the Higgs field, ($} = (0, v) J , the SU(2) L - 
, doublet neutrinos vn, receive a non-zero mass matrix of the form (M„ L )ij = XijV 2 /A. 

^? . There exist three tree-level realizations of this effective dimension-5 operator using only renormalizable in- 

teractions. The simplest one, widely known as the type-I seesaw Q, requires the extension of the SM by a number 
hr singlet fermions, which are usually taken to be heavy right-handed (RH) neutrinos v a R (with a = 1, 2, ... , Hr). 



After integrating out the heavy singlet neutrinos v a R, we obtain Weinberg's effective dimension-5 operator mentioned 
above. Specifically, the neutrino masses and mixings may be deduced from the Yukawa Lagrangian: 

- £y = {Vu)iaLi^V aR + ^V < ^ R {M Ii ) al 3Vf3R + H.C. , (1) 

which in turn implies the following (3 + Ur) x (3 + tir) complex symmetric matrix in the flavor basis {(^l) i v oirY- 

*• - ( k Z ) ■ 

where Mr> — y v v is the Dirac mass matrix for the neutrinos and Mr is the (B — L)-breaking Majorana mass matrix 
of the RH neutrinos. In the usual seesaw approximation (Mr})i a /(MR) a p <C 1, the physical light neutrino masses are 
given by 

M VL ~ - M D M^Ml . (3) 

Evidently, the seesaw scale Mr is correlated with the magnitude of the Yukawa coupling y v . For the experimentally 
observed sub-eV light neutrinos [![, Mr is usually constrained to be much larger than the electroweak scale mw, 
unless y v is very small (5s 10 -6 ) [cf. (jSJ)] or there are cancellations in the mass matrix structure given by © due to 
specific flavor symmetric patterns of y v and Mr Q. 

An interesting realization for the seesaw scale to be in the TeV range is the so-called inverse seesaw model Q, 
where in addition to RH neutrinos {y a R\, another set of SM singlet fermions {S p l} are introduced. Here we consider 
a symmetric extension of the SM with three pairs of singlet neutrinos (i.e. a, p = 1, 2, 3). One of the salient features 
of the inverse seesaw mechanism is that the small neutrino masses are generated by a small lepton-number breaking 
3x3 mass matrix fis- In contrast to the type-I seesaw case, the mass matrix of the observed left-handed neutrinos 
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(f eijL ( )T )r vanishes in the limit of lepton-number conservation: fig — > 0. More explicitly, the neutrino Yukawa sector 
of a general inverse seesaw model is described by the Lagrangian 



-C Y = (y v )ia Li$v aR + S pL (M N ) pa v aR + -i>2 R {nR) a pvpR + 7; s pl{vs) p \Sxl + H.c. , (4) 



which gives rise to the following 9x9 neutrino mass matrix in the basis {(vil) > v aRl (S p l) c }- 



( M D \ 

M v = Ml m Ml . (5) 
\ M N ii S J 



Note that we have not included in Cy the dimension-4 lepton-number breaking term L^Sj^ which appears, for 
instance, in linear seesaw models B, since the resulting neutrino mass matrix in presence of this term can always 
be rotated to the form given in ([5]) !9j. Observe that the standard inverse seesaw model discussed originally in Q is 
recovered, once we set the RH neutrino Majorana mass hr — in ([5]). The smallness of the 3x3 symmetric matrices 
fJ-R,S are 'technically natural' in the 't Hooft sense. In other words, in the limit of hr,s ~~ ^ 0, lepton number symmetry 
is restored and the light neutrinos V\,i, 3L are massless to all orders in perturbation, whereas {(2^1,2,3)0 ., (5*1, 2, 3)^} 
form three heavy singlet Dirac neutrinos of approximate mass (Mn) 1 2 3 . 

In this paper, we analyze another interesting realization of inverse seesaw models, where \ir ^ 0, but jj,g = 0. 
In this case, the light neutrinos are massless at the tree level, but acquire a small mass at the one-loop level. We 
show that this one-loop induced light neutrino mass is directly proportional to the Majorana mass matrix /ir and 
results from well-known SM radiative corrections that involve the Z- and Higgs bosons pj|. We call this scenario, the 
Minimal Radiative Inverse Seesaw Model. This is a rather economical scenario as it does not require the existence of 
other non-standard scalar or gauge fields or other fermionic matter beyond the singlet neutrinos {v a R, S p l}- Finally, 
we also analyze the general case, where both [ir and (is are non-vanishing and present numerical estimates regarding 
the typical relative size of the two Majorana mass matrices /ir and /xs, while maintaining agreement with neutrino 
oscillation data. 

This paper is organized as follows: In Section II, we discuss the tree-level neutrino mass matrix in inverse seesaw 
models. In Section III, we present the one-loop neutrino masses and mixings in the Minimal Radiative Inverse Seesaw 
Model. In Section IV, we analyze the general inverse seesaw model at one-loop level. Finally, our conclusions are 
given in Section V. 



To obtain the physical neutrino masses, we cast the mass matrix in ([5]) into a form similar to the type-I seesaw 
case given in ([2} by defining a 3 x 6 mass matrix A4d = (Md,0) and a 6 x 6 mass matrix Ais in the weak basis 



We may now perform a block diagonalization of the 9x9 neutrino mass matrix ([5]) by a unitary transformation V to 
cast it into the form: 



II. TREE-LEVEL NEUTRINO MASSES 



{K2,3)*,0Sl,2,3)£}: 




(6) 




(7) 



where the unitary matrix V has an exact representation in terms of a 3 x 6 arbitrary matrix £ |ll| 




(8) 



giving rise to the exact analytic expressions for the block mass eigenmatrices [l2j 1 : 

m ur = (i 6 +^o- 1/2 (m s +m t d c+^m d )(i 6 +c t c)- 1/2 , (9) 
m vl = -(i 3 + cc r r 1/2 (M D c T + (m t d - (Ms( T )(i 3 + cc T r 1/2 = -cm„ r c t . m 

The block diagonalization condition implies that the (1,2) block element of V T A4„V in ([7|) vanishes, i.e., 

M D - (M s - (MlC* = 3 x6, (11) 

which can be solved for ( in terms of Md and Mg 1 . Moreover, M.'g 1 is obtained by inverting the mass matrix given 
in ©: 

M -i = ( (hr - mI^MnY 1 - {hr - mI^Mn)' 1 mJ^ 1 \ (12) 

Note that Mg 1 given by (|12[) is symmetric, since it is easy to show that for symmetric invertible matrices /t_R,s, 

(^s-MN^Miy 1 Mnh^Y = (fiR -MIh^Mn)' 1 Mjffig 1 . (13) 

The physical neutrino mass matrices given by (|9]) and (|10[) can be expanded in a converging Taylor series in £, 
provided its norm is much smaller than unity, i.e., ||£|| = y/Tr(0() <C 1, which is naturally satisfied within the generic 
seesaw framework. To leading order in (jllj) implies £ = Md-MI^ 1 , and hence the light neutrino mass matrix 
given by (|10[) simplifies to 

M VL = -M D M s l M J D + 0(\U 2 ) 

= -M D {M s l ) URU cM T D + 0{\\C\\ 2 ). (14) 

This is analogous to the type-I seesaw formula (JU), which is shown diagrammatically in Figured] Notice that it is 
only the vrv'r component of Mg 1 that contributes to the light neutrino mass matrix at the tree level. 

(*> (*> 
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FIG. 1: The tree- level diagram for the light neutrino mass in the Inverse Seesaw Model. 

In the limit ||/iR,s|| ^ ||-Mjv||; the (1, 1) block of .M^ 1 given by (|12j) can be expanded in powers of /is. We then 
obtain from (fl4|) the tree- level light neutrino mass matrix 2 : 

MX CC = MdM^Hs {Ml)' 1 Mj, + MdM^hs (Ml)' 1 UrM^hs (M^)' 1 M t d + 0( f i 3 s ) , (15) 

which identically vanishes in the limit /is — > 0. This can also be deduced from the exact expression in (|10[) . since in 
this limit, C,ftA T D = which implies C, = MdM.^ 1 from (fTT|) . and hence, M VL = in (| 10|) . Therefore, in this limit, 
the full neutrino mass matrix M v in ([5]) has rank 6 and the light neutrinos remain massless at the tree level, even if 



1 Some approximate seesaw expressions are also studied in Il3tl ■ 

2 Notice that 1 1 /ir g \ | could be bigger than 1 1 M r> \ \ in this inverse seesaw approximation. 
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FIG. 2: One-loop diagrams pertinent to the generation of light neutrino masses in the Minimal Radiative Inverse Seesaw Model. 



III. RADIATIVE NEUTRINO MASSES AND MIXINGS AT THE ONE-LOOP LEVEL 



At the one-loop level, the light neutrino mass matrix M VL acquires a small radiative mass from the self-energy 
diagrams shown in Figure [5] We have drawn the figures in the flavor basis to explicitly show the lepton number 
violating mass insertions. Note here that the scalar propagator $ includes both the neutral Higgs boson (H) and 
the neutral Goldstone boson [G ) contributions. In the on-shell renormalization scheme and in the Feynman gauge, 
these are the only one- loop diagrams contributing to the neutrino mass, as in the type-I seesaw case [lol [T3 . [l5j . 
In a general i?^-gauge parameterized by the gauge-fixing parameter £z, the £z dependence of the self-energy for the 
Z-loop diagram is exactly canceled by the G°-loop diagram, and the final result is independent of £z [HI- 

The one-loop induced light neutrino masses are obtained from the momentum-dependent self-energy as follows: 



M l-loopp fl 



(16) 



where Pr = (1 + 7 5 )/2 is the right-chirality projection operator. Here Yi(jf) is the sum of the self-energies of the 
diagrams shown in Figure [H and is given by the following gauge-invariant Feynman amplitude: 



p R mMPR = 



9 



4 cos 2 



-Mi 



9 ™H 
4m^ 



d d k 



Mi 



d d k 



Pr 



■Mi 



(27T) 



(f - t) 



(p — k) 2 + ie 



k) 2 l 6 -M 2 s +ie 
(f-f)U 



_(fc 2 



Mi 



(p-fc) 2 le 



M% 



1 



ie) 2 



PrM t d 



(j) — k) 2 + ie 



1 v Pr 
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(17) 



M\, 



D > 



where g is the weak gauge coupling, W is the weak mixing angle, and d = 4 — e is the dimensionality of space-time in 
the dimensional regularization scheme. Notice that the integrands in (|T7|) vanish in the infra-red (IR) limit fc M — > 0, 
because of the tree-level vanishing condition on light neutrino masses. However, the loop momentum k is integrated 
over all possible values, and so we get a non-zero contribution to the light-neutrino mass matrix at the one-loop level. 
A straightforward evaluation of the ultra-violet (UV) finite Feynman integrals appearing in (|17l) leads to the following 
neutrino mass matrix in (1161): 



M 1_loop 



aw 



167rra^ 



M D M S 



M% 



In 



Ml 



•h 



3m| 



Ml 
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(18) 
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FIG. 3: The analytic dependence of the one-loop function }{xm) [cf. (|20f) ] and a;jv/(a;jv) on the heavy neutrino mass scale xm- 



where a w = g 2 /An is the weak fine structure constant. In the limit ||//r,s|| <C ||Mjv|| and assuming that Mjv = mpfl^, 
for simplicity, we derive the simpler expression 



M l-loop 

VL 



where the one- loop function /(xjv) is given by 



2~ M D fi R M D 

m w 



(19) 



I{ x n) 



aw 
16tt 



x H 



XN - X H 



In 



3a; z 



XN ~ X Z 



hi 



xz 



(20) 



with a; at = m N /m' w , xh = m H /m w , xz = rn 2 z /m w = 1/ cos 2 6*„,. From P^l) . we clearly see that at the one- loop 
level, the light neutrino mass depends linearly on the lcpton number breaking term /in (of dimension-3), and hence it 
does not vanish in the limit /is — 0> unlike at the tree level. 

The one-loop function /(xjv) given by (|2T)|) is plotted in Figure^ for xh — 1-88 (corresponding to run — 125 GeV), 
xz = 1.30 and aw = 0.033. Note that in the limit xn 3> 1 (as in the canonical type-I seesaw), the one-loop function 
stated in (f2"U)) becomes 



aw xh + 3xz 



ln(x]v) 



167T XN 

Hence the light neutrino masses given by (|19[) are suppressed by l/m N : 



M 



l-loop ^ a W M D y, R Mj) m 2 H + 3m| 



16tt 



In 



''N 



'■w 



N 



'W 



On the other hand, it is interesting to note that in the limit xjv *C 1, the one- loop function takes on the form: 



J{xn) 



aw 

167T 



In 



V x 



This leads to a light neutrino mass which is only logarithmically dependent on the heavy neutrino mass: 

vfi-loop _ M D fi R Ml ( m H m% 

~ 16tt 



I 4 



(21) 



(22) 



(23) 



(24) 
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It is important to stress here that this ln(mjv) dependence of the light neutrino mass is a unique feature of the Minimal 
Radiative Inverse Seesaw Model we have been studying here. 

Finally, we wish to offer some comments on the validity of our one-loop analysis presented here. In a general study, 
one should also consider the loop corrections to the other mass terms appearing in the full inverse seesaw mass matrix 
in ([3]), as well as wave-function and mixing renormalizations [Til. Il6j. However, in the MS renormalization scheme, 
the UV divergences of the Dirac mass in the v^-vr sector can be absorbed in the renormalized Mo- On the other 
hand, all singlet Majorana and Dirac mass matrices in the vr-Sl sector do not get any UV-divergent contributions, 
at least at the one-loop level. 



IV. THE GENERAL INVERSE SEESAW MODEL 



Let us now consider the general inverse seesaw scenario, where both the Majorana mass matrices /is and /j.r do 
not vanish, i.e. ^5 7^ 0, \ir 7^ 0. To leading order in [Ir.s, the light neutrino mass matrix is given by 

M VL = ^M D Ls+XNf(x N )p R )Ml + 0(i4,iJ, 2 R ) , (25) 

where the one- loop function /(xn) is defined in (|20|) . 

In order to get an order of magnitude estimate of the relative size of [1r with respect to /is required to fit the 
neutrino oscillation data, let us choose a basis in which the charged lepton mass matrix is diagonal. Then we may 
write the light neutrino mass matrix as 

M VL = U T M^U , (26) 

where U is the PMNS mixing matrix 3 given in terms of the three mixing angles, one Dirac and two Majorana CP 
phases [l|: 



C12C13 S12C13 s 13 e l " 



U = -S12C23 - c 12 s 23 s 13 e i& c 12 c 23 - s 12 s 23 s 13 e u s 23 c 13 | x diag(e iQl / 2 , e 1 ^' 2 , 1) , (27) 
V S12S23 - c 12 c 23 s 13 e iS -c 12 s 23 - s 12 c 23 s 13 e lS c 23 c Vc 

with Cij = cos 6ij , Sij = sin 0,j . For a normal hierarchy, one has 



= diag(m 1 , y'Am^, . (28) 

Using the central values of a recent global-fit analysis of the neutrino oscillation parameters [l8j |: 

Am s 2 ol = 7.62 x 1CT 5 eV 2 , Am 2 tm = 2.55 x 1CT 3 eV 2 , 12 = 34.4°, 6 23 = 40.8°, 13 = 9.0°, 6 = 0.8tt , (29) 

and assuming mi — 0,ai = a 2 — 0, we obtain from (126[) : 

/ 0.31-0.12i -0.09 + 0.32i -0.72 + 0.37i \ 
M VL = (10^ 2 eV) x -0.09 + 0.32i 2.53 + OMi 2.19 + O.Oli . (30) 
\-0.72 + 0.37i 2.19 + 0.01i 3.07 - 0.03i / 

For illustration, let us absorb the flavor structure of Mb in (|25|) into the matrices /iR.s, by rewriting the light- 
neutrino mass matrix as 



A_d \Us + x N f(x N )fiR^j = XdUcS , (31) 



3 For simplicity, we have not considered here the leptonic non-unitarity effects which are generally of order 1 1£| | in inverse seesaw models |17| . 
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where X]j = m 2 D /mj^ and m-D max = max (l(MD)ij|) is the largest element of (in absolute- value terms). In 
addition, we have re-defined ^s,R as 

fl s = Kfi s K T , fx R = K f i R K T , (32) 

where K = MD/mD max is in general a dimensionless 3x3 complex matrix. With these definitions, we can now make a 
one-to-one mapping between the elements of M VL and p.R,s- This is shown in Figure U where we have plotted against 




FIG. 4: The magnitudes of the largest element of us (dashed line) and (solid lines for various values of xn) as a function 
of the ratio \d ■ 

the ratio Ad: (i) the magnitude of the largest entry in fts (denoted as jl'g ax ) for JIr — (the dashed line); (ii) the 
magnitude of the largest entry in p,R (denoted as /2™ ax ) for jls — and for xn — 0.01, 1, 100 (the solid lines). We 
note that the mass scale of /i"g ax is roughly 2-4 orders of magnitude larger than that of //™ ax . Thus, if /jr ^ 0, a 
much milder hierarchy can be realized between the lepton-number breaking scale /i^ ax and the electroweak scale mw 
in the Minimal Radiative Inverse Seesaw Model. 

In the general inverse seesaw scenario where (Ir^s Oi the mass matrix p, c f{ defined by (|31|) : 

Pes = fis + x N f{xN)p-R, (33) 

is fixed by the neutrino oscillation data for a given value of Ad. The dependence of £t™g* on Ad is exactly the same 
as that of /I™ ax (with Jjlr = 0), as shown by the dashed line in Figured] However, the relative size between hr and 
/is depends on their relative sign, for a fixed given ratio xn- For instance, when both /I™ ax ,/i^ ax > 0, the constant 
contours of /iSf* (in keV) are shown in Figure [5] for = 1. For /2™ ax < and /lg ax > 0, the corresponding contours 
are shown in Figure [H where we have the cancellation regions (/i"Jf x — > 0) for certain combinations of (/2g lax , P'r^)- 

We conclude this section by commenting that our numerical results only depend on the largest element of the light 
neutrino mass matrix which is always of the order of y/ Am^ tm , irrespective of the neutrino mass hierarchy. Hence 
our order of magnitude estimates will be valid for an inverted hierarchical light neutrino mass spectrum as well. 



V. CONCLUSIONS 



We have presented a minimal radiative mechanism for generating light neutrino masses in inverse seesaw models. 
The radiative neutrino masses arise at the one-loop level from known SM electroweak quantum effects involving 
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FIG. 6: Contours of constant /2™ ft ax (in keV) for x N = 1 and for fig** < 0, /i£ ax > 0. 



the Z and Higgs bosons. Unlike in other radiative inverse seesaw mechanisms existing in the literature [19(, the 
implementation of our radiative mechanism docs not require any extra non-standard fields other than usual SM 
singlet neutrinos. 

In the Minimal Radiative Inverse Seesaw Model, where the dimension-3 lepton-number breaking mass matrix /i^ of 
the right-handed neutrinos is non-zero, we have found that the light neutrino masses generated at one-loop level are 
UV- finite and are directly proportional to [ir. We also showed that the 3-by-3 Majorana mass matrix /ir could be 2 — 4 
orders of magnitude larger than the other 3-by-3 Majorana mass matrix /is present in the standard inverse seesaw 
models. Hence, this could alleviate the hierarchy between the lepton number breaking scale and the electroweak scale 
in these models. In a supersymmetric version of this mechanism (e.g., (20j). one might expect to ameliorate this 
hierarchy even further, in addition to having a scalar dark matter candidate in the form of the lightest sneutrino (2lj . 
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We hope to return to these issues in a future communication. 
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